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Abstract 

We report experimental and numerical analysis of expontentially shaped long Josephson junctions 
with lateral current injection. Quasi-linear flux flow branches are observed in the current-voltage 
characteristic of the junctions in the absence of magnetic field. A strongly asymmetric response 
to an applied magnetic field is also exhibited by the junctions. Experimental data are found 
in agreement with numerical predictions and demonstrate the existence of a geometry-induced 
potential experienced by the flux quanta in nonuniform width junctions. 

PACS numbers: 74.50. +r 
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I. INTRODUCTION 



In recent years, there have been some theoretical studies concerning the possibihty to 
influence the flux motion in long Josephson junctions by mean of geometry- induced or field- 
induced potentials. The most known example is the annular iunction embedded in a spatially 
,o„. „.g„et,c field. A. ,t . U„ow„H, ,„ tMs case a cosinusoidal potential . expe- 
rienced by a flux quantum trapped in the junction when a spatially homogeneous magnetic 
field is applied parallel to the junction barrier. The origin of the potential is caused by the 
spatial variation of the radial component of the magnetic field that in this special geometry 
becomes sinusoidal. A field-induced sawtooth-like potential has been recently consideredi for 
experimental demonstration! of ratchet effect in annular junctions. Currently, modifications 
of the annular geometry, as the heart-shaped geometry! are used to achieve a field-induced 
double-well potential for demonstration of fluxon quantum-bit and macroscopic quantum 
coherence phenomena. 

To achieve potentials without the help of a magnetic field, the case of a nonuniform 
junction width has been theoretically addressed in recent years for linearli'lil as well as 
for annular! geometries. Theoretically, a geometry-induced potential related to the spatial 
variation of the junction width is expected. This potential corresponds to a force acting 
on the fluxons in the direction of the shrinking width. Recently,^ ordinary Josephson flux- 
flow oscillatorjnilli have been modified adding to the classical overlap geometrylli unbiased 
pointed tails. This is expected to enhance the annichilation of the fluxons at the edges of 
the oscillator, with consequent reduction of fine structures in the ordinary velocity-matching 
step0'00. 

In this paper we experimentally address the existence of the geometrical force in nonuni- 
form width junctions. To do this, we consider an exponentially shaped overlap junction with 
lateral current injection. The lateral current injection acts as a flux quanta generator, also 
in the absence of a magnetic field, while the unbiased shaped region should act as an accel- 
erating region for both fluxons or antifluxons injected at one edge of the junction. If really 
present, the geometric force should allow to achieve a quite regular unidirectional flux flow 
motion in the junction, without the help of a magnetic field. This dynamical regime should 
be accounted for a branch in the current- voltage curve of the shaped junction. The demon- 
stration of such branches in the current-voltage curves of shaped junctions we fabricated 
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suggests that such a geometrical force is really experencied by the flux quanta. 

The paper is organized as follows. In Sec. II we specialize the general mo dei'B for a long 
junction with nonuniform width to our exponentially shaped junction with lateral current 
injection. In Sec. Ill the experimental results both in the absence of magnetic fleld and 
in the presence of magnetic fleld are presented and discussed with the help of numerical 
simulations. Finally, main results are summarized in the Conclusions. 



II. THEORY 

Under some simpliflcative hypotheses!, the model for a long overlap junction with nonuni- 
form width W{x) = fi{x) — f2{x) was found as 



sin 6 + ac 



with 



T(x) 



W'{x) W'jx) 



W{x) 



r(x). 



(1) 



(2) 



In Eq. (|l]) is the Josephson phase, a is the dissipation parameter, rj^ and rjy are the 
normalized magnetic flelds in the x and y directions, respectively. Space is normalized 
to the Josephson penetration length Aj and time to the inverse of the plasma frequency 
ujj = c/Xj, with c the velocity of electromagnetic waves in the junction. For the geometry 
we report here (see Fig. 0) the total physical length of the junction is L = Lq + Ls + Ll ^ Aj, 
while the width is chosen as 



W{x) 



Wo <x<Lo 
Wq exp 

Wl 



Ls + Lq < X < L, 



Lq < X < Ls + Lq 



(3) 



The bias current I can be fed into the left or into the right edge, as shown in Fig. ^ The 
bias term T{x) in Eq. becomes 



Lo < X < L, 



(4) 
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FIG. 1: Exponentially shaped junction with lateral current injection. The current can be fed into 
the left [(a)] or into the right [(b)] edge. 

for left edge current injection, and 

, , |o 0<x<Ls + Lo, 
1b{x) = < ^ (5) 

Ijolffel^^^- Ls + L,<x<L, 

for right edge current injection. The lengths Lq, Ll, Wo,and Wl are chosen to be shorter 
than Aj and such that LlWl — LqWq, to have equal biased areas. 
Hence, the model for our exponentially shaped junction becomes 

(f^xx - (t>tt = sin + a(t)t + Mx -r]\- ^a,b{x), (6a) 
ifM = V, (6b) 
^x{l) = V, (6c) 

where / = L/Xj, A = Xj\n{Wo/WL)/ Ls, and r] accounts for an external magnetic field 
applied in the y direction. We remark that the chosen geometry reduces to the exponentially 
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shaped in-hne geometryou in the case of left current injection with injection length Lq much 
lower than Aj. In such a limit, the model becomes 

4>xx - 4>u = sin + a(f)t + X4>x - tiX, (7a) 
¥^.■(0) = V-lh (7b) 
Vx{l) = V- (7c) 

As first noted in Ref. a force that drags fiuxons or antifiuxons in the direction of the 
narrowing width is expected for our geometry. In fact, in the absence of magnetic field 
[r] = 0), for a soliton 



(f){x) 



4 arctan 



exp 



X — ut 



a 



in the unbiased region Yia.b{x) = 0], the following equation of motion can be found from 
Eqs. (|3) in the framework of the perturbative approachlll 



(1 - uY'^'^ = -"^^ + -F^- (8) 

This indicates that, if present, both a fiuxon or an antifiuxon will experience a force propor- 
tional to the shaping parameter A and will be accelerated in the narrowing width direction, 
i.e. from the left to the right in our geometry. From Eq. (|^) the stationary velocity of the 
motion will be m = A/a for A/a < 1 or the limit velocity m = 1 for A/a > 1. 



III. NUMERICAL AND EXPERIMENTAL RESULTS 

We realized Nb/Al203/Nb junctions with the geometry shown in Fig. |I[ The physical 
dimensions of the junctions were L = Lq + Ls + = (10 + 560 + 40) /zm, Wq = 40 /im, 
Wl = 10 fim. For the two junctions we report here the normalized lengths were I ~ 20, 
and / ~ 17, with shaping parameters A ~ 0.07 and A ~ 0.08, respectively. In the following 
the behavior of the junctions in the absence of magnetic field as well as the response to a 
magnetic field applied along the y direction is discussed. 
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FIG. 2: Current-voltage curve of a junction with left current injection (a) or with right current 
injection (b) exhibiting quasi-linear flux flow branches. In (c) the flux flow branch for left current 
injection is compared with the one achieved for right current injection. 

A. Behavior in the absence of magnetic field 

The current- volt age curve for the junction with / 20 is reported in Fig. |^. No external 
magnetic field is applied to the junction. Panel (a) of this figure refers to the case of current 
injected into the left edge, while the panel (b) refers to the case of current injected into the 
right edge. In both cases an almost linear branch starting from a certain critical current 
and extending for a given current range is observed, but some qualitative differences exist 
between the two cases. As better seen in Fig. 0(c), when current is injected into the left edge 
the branch starts to a lower critical current and is more regular than the branch achieved with 
current injected to the right edge. The observed behavior seems to be consistent with the 
idea that a geometrical force is effectively experienced by the fluxons in this geometry. When 
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current of positive (negative) polarity is injected into the left edge, antifluxons (fluxons) will 
enucleate at this edge after a critical current value has been reached. From Eq. (|^) these 
antifluxons are expected to be accelerated toward the other edge because of the geometrical 
force, so realizing a regular unidirectional flux motion. Conversely, if the current is injected 
into the right edge, the antifluxons enucleated a this edge should overcome a force to travel 
toward the left edge. This should result in an irregular or chaotic flux motion. Moreover, 
in the left current injection the geometrical force helps to start an antifluxon (or fluxon) 
motion, corresponding to a voltage in the I-V curve beyond a critical current value. In the 
right current injection such a force opposes the starting of the flux motion, so resulting in a 
larger critical current. 

We should remark that the quasi-linear branches reported in Fig. ^ could remind the 
Displaced Linear Slope branches sometimes reported for rectangular in-line or overlap 
geometries0'0'@ffl. However here the branches are rather regular, quite noiseless, and 
are obtained in the absence of magnetic field. 

In the following we will focus on the left current injection. In Fig. ^(a) same data of 
Fig. |^(a) are replotted on a larger scale. As better seen in the inset, the flux flow branch 
exhibits a series of small steps spaced of a voltage AV ~ 15 fiV. The curve reminds the 
numerically predicted curve for an exponential shaped asymmetric in-line junction!. The 
observed voltage spacing is consistent with the spacing of cavity mode resonances expected 
from physical dimensions of our junction, AV ^qc/L. This is to be expected due to the 
open boundary condition at the edges. The antifluxons in the chain moving toward the right 
edge will be reflected as fluxons. If the dissipation a is not too large, the antifluxons can have 
enough energy to travel toward the left edge as fluxons after the reflection, despite of the 
geometrical force opposing the motion. This mechanism can excite cavity mode resonances. 
However, for quite large dissipation one expect that the reflected motion would be more 
and more damped, and the excitation of cavity modes should be consequently damped. In 
Fig. ^(b) the I-V curve of the junction with / ^ 17 is plotted for three different temperatures, 
corresponding to three different a values. As expected, the small steps accounting for the 
cavity mode resonances are more and more damped as the dissipation (temperature) is 
increased. 

To gain further insight in the flux motion dynamics in the absence of magnetic field, we 
integrated the model Eqs. (|^) with rj = and forcing term 7a (a^) defined in Eq. (^). In 
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FIG. 3: (a) The current-voltage curve of Fig. |2|(a) is replotted on a larger scale. Small steps in 
the flux-flow branch are shown in the inset, (b) Flox flow branch of a junction with A w 0.08 at 
different temperatures. 

Fig. H(a) we show the calculated current voltage curve for a junction with left edge current 
injection and of uniform width [A = in Eqs. @]. This is equivalent to the asymmetric 
in-line rectangular geometrylli. As seen in the snapshot showing the instantaneous voltage 
distribution in the junction, antifluxons are created at the biased edge. However, due to the 
absence of a force in the unbiased region, the flux motion is not very regular. Here it is only 
the repulsion between flux quanta that tends to drive the flux toward the right edge. The 
resulting motion is quite noisy, as well as the calculated ac voltage at the right edge. 

In Fig. ^(b) the case of an exponentially shaped width (A = 0.07) is considered. In the 
simulation we used parameters similar to the ones estimated for the experimental curve in 
Fig. 0(b). As it seen in the snapshot, now the presence of a force in the unbiased region 
makes the flux motion toward the right edge more regular, as well as the voltage signal at 
the right edge. In the snapshot four moving antifiuxon are counted as present in the mean 



8 




FIG. 4: Calculated current-voltage curve for an unmatched rectangular [(a)], for an unmatched 
exponentially shaped [(b)], and for a matched exponentially shaped [(c)] junction with lateral 
current injection. The snapshots shows the instantaneous voltage profiles and the voltage signal 
at the right edge of the junctions. 

in the junction for the chosen bias point. By increasing the bias current, more and more 
antifiuxons can be injected in the junction, and the chain becomes more and more dense. 
Correspondly, the voltage signal at the right edge becomes less impulsive and approaches a 
sinusoidal form. When the junction is completelly filled with antifiuxons, a transition to a 
new dynamical regime, similar to a laminar phase fiowi is achieved. In the current-voltage 
curve, this transition corresponds to the switch from the fiux fiow branch toward the other 
resistive branch, beyond a critical current value. As in the experimental curve, the small 
steps in the calculated fiux fiow branch are spaced of tt/Z, the spacing (in normalized units) 
expected from cavity mode resonances excited from refiecting boundary conditions at the 
edges. 

As said above, the resonances can be damped by increasing the dissipation a. However, 
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another mean to achieve this is to match the impedance of the fluxon chain with a load z 
at the right edge. The matching is found possibleB'i for this geometry, conversely to the 
rectangular asymmetric in-line geometry. For our exponentially shaped junction the wave 
impedance is just the stationary velocity we have found above, —(pt/'Px = u = \/ a. This 
should be matched with a load of impedance z = —(f)t{l)/(f)x{l). Numerical results for the 
case of a matched load are shown in Fig. ^(c). As a result of the absence of reflections, the 
small steps typical of the unmatched case of Fig. @(b) are now absent and the flux chain 
exhibits a very regular motion, as well as the voltage signal at the matched edge. This 
characteristic makes the exponentially shaped junction interesting as zero-magnetic field 
flux flow oscillator. 

In Fig. ^ there are shown numerical results for a junction with A/a > 1. For the unmatched 
case shown in Fig. ^(a) the used parameters can account for the current-voltage curve at 
T=4.2 K in Fig. ^b). In both cases A/a > 1. As said above, in such a case the antifiuxons 
are accelerated toward the asymptotic velocity -u = 1, so the matching condition is now 
z = 1. In the snapshots shown in Fig. ^(a) and (b) we show also the instantaneous magnetic 
field. Antifiuxons in the chain are accelerated and, due to their relativistic nature, they are 
Lorentz-contracted as the asymptotic velocity is approached. As for the case A/a < 1, the 
matched junction exhibits a smooth fiux-fiow branch and the voltage signal at the right edge 
is very regular. 

B. Behavior in the presence of magnetic field 

Due to the lateral current injection and due to the existence of a preferred direction of 
motion, our exponentially shaped junction is expected to show a behavior in magnetic field 
even more strongly asymmetric with respect to the asymmetric rectangular in-line geometry. 
The calculated critical current as a function of the magnetic field is shown in Fig. |^(a) for a 
junction with / = 20 and A = 0.07. The pattern was obtained integrating Eqs. (^) with 1] 0. 
A quite abrupt decrease of the critical current around rj = 2 is recovered. In normalized 
units, for this value of magnetic field a fluxon is enucleated in the junction. The almost 
complete absence of secondary lobes in the pattern means the almost complete absence of 
trapped flux in this kind of junction, an indication that the geometrical force helps to move 
fluxons just after their enucleation. 
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FIG. 5: Calculated current-voltage curve for an unmatched [(a)] and a matched [(b)] exponentially 
shaped junction with A = 0.12. The snapshots shows the instantaneous voltage and magnetic field 
profiles, and the voltage signal at the right edge of the junctions. 

Figure P(b) shows a global representation of calculated resonant steps achieved for nor- 
malized magnetic field values slightly larger than the critical value rj = 2. Such a plot is 
obtained superimposing the curves corresponding to different values of magnetic field. Two 
families of steps with two characteristic voltage spacing are recovered. The lower voltage 
spacing is the one expected for Fiske mode resonances, A ~ n/l, the larger one is about 
27r// = 2 A. The family of steps with larger voltage spacing appears in the same current 
range where the flux flow branch is achieved in the absence of magnetic fleld. From numeri- 
cal simulation it is seen that these steps with larger voltage spacing consists of cavity mode 
resonances excited by a fluxon chain and an antifluxon chain moving in opposite directions. 
In fact, for positive bias current an antifluxon chain moving toward the rigth edge is gener- 
ated, while a positive magnetic fleld generate fluxons. These fluxons are pulled toward the 
left edge by the Lorentz force associated to the positive bias current a this edge. Antifluxons 
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FIG. 6: (a) Calculated magnetic field pattern for a junction with A = 0.07. (b) Calculated steps 
induced in the current voltage characteristic by a magnetic field varied in the range 2.5 < rj < 3.3. 
In the inset the steps at for a given field value are shown. 

travelling toward the right generate a voltage with positive polarity that adds to the voltage 
of the fiuxons traveling toward left. This accounts for the larger voltage spacing observed 
for the steps of this family. 

For negative bias currents, the fiuxons generated by the positive magnetic field are pushed 
toward the right edge by the Lorentz force associated to the negative bias current. Moreover, 
also the fiuxon chain injected at the left edge by the bias current moves toward the right, 
under the action of the geometrical force. The result is that an unidirectional motion toward 
the right edge tends to be achieved, with consequent damping of resonant steps accounting 
for cavity mode resonances. As seen in the inset of Fig. ^(b), resonant steps are in fact 
virtually absent for negative current values. This peculiarity is observed also for larger 
magnetic values. For negative magnetic fields values, the current-voltage curve show steps 
for negative bias current values, and no steps for positive bias values. In other words, the 
current-voltage curve reflects with respect to the origin when the sign of the magnetic fleld 
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FIG. 7: Measured magnetic field pattern (a) and current steps (b) induced in a junction with 
A ~ 0.07 by a magnetic field slightly larger than the critical field. In the inset the current voltage 
curve at B=2.4 G is shown. 

is inverted. 

Predicted strongly asymmetric magnetic field behavior summarized in Fig. |^ is fully 
recovered in the experiment, as shown in Fig. 0. Experimental data refer to the junction 
with I ^ 20. In Fig. ^(b) the steps are achieved for magnetic field values slightly larger than 
Be = 1.75 G, the critical field of the junction. In Fig. |^(a) we reported the modification of 
the fiux fiow branch induced by a low magnetic field, i.e., lower than B^. An asymmetric 
tunning of the branch is recovered . This can be easily understood as follows. We are using 
strongly left-edge-peaked current injection. In this case our junction can be also described 
as a shaped aysmmetric in-line junction, i.e., with model Eqs. (^. Looking at the vortex 
generator term ^^^(O) = 77—7/, it is easily understood that a positive magnetic field cooperates 
with a negative bias current to inject fiuxons in the junctions, while such a positive magnetic 
field opposes a positive current injecting antifiuxons in to the junction. The result is that, 



13 




FIG. 8: (a) Modification of the flux flow branch induced by a magnetic field lower than the critical 
field, (b) Steps recorded for magnetic fields quite larger than the critical field of the junction. A 
smooth single step is achieved for B > 8 G. In the inset we show the voltage of the junction 
polarized at fixed current on the step as a function of the magnetic field. 

for a given absolute value of the bias current, more solitons are present at negative polarities 
than at positive polarities, resulting in a voltage at negative polarity larger than the voltage 
at positive polarity, as it is in fact observed in the experimental curves of Fig. |^(a). 

Figure 1Kb) shows the steps recorded in the junction for magnetic fields quite larger than 
Be- The fine structures accounting for Fiske modes are completely smeared out at i? ^ 9 
G and a smooth step, similar to the velocity-matching step000, with asymptotic voltage 
strictly proportional to B is achieved. The proportionality between the voltage and the 
magnetic field is shown in the inset for a given biasing current. The velocity-matching step 
is found to exist beyond the current- voltage range of existence of the linear flux flow branch 
in zero field, i.e., beyond V =600 /iV (and up to 1500 fiV) and / = 9.5 mA. This is consistent 
with the observation0 that the velocity matching step originates from a quasi-linear waves 
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regime in the junction. In our case such a quasi-hnear background is the laminar phase 
achieved beyond the current range typical of the zero-field fiux-fiow branch. 



IV. CONCLUSIONS 

Summarizing, we have experimentally investigated the occurence of dynamical states in 
exponentially shaped overlap junctions with lateral current injection. In zero magnetic field, 
the lateral current injection acts as a fiuxon or antifiuxon chain generator and these chains 
can be accelerated by the geometrical force originating from the nonuniform width of the 
junction. The result is that a quite regular fiux-fiow motion, corresponding to a quasi-linear 
branch in the I-V curve of the junction, can be achieved in this kind of junction without 
the help of an external magnetic field. Moreover, numerical simulations show that this kind 
of motion can be precisely matched to a load, a peculiarity that makes the zero field fiux 
fiow branch interesting for zero-field fiux fiow oscillators. In the presence of a magnetic 
field a rather asymmetric behavior is exhibited by the junction. For low magnetic fields an 
asymmetric tuning of the fiux-fiow branch is observed, for moderated magnetic field ordinary 
Fiske modes steps are mixed with nonlinear cavity modes steps, and for quite large magnetic 
fields a single step similar to the velocity-matching step known for uniform width geometries 
is recovered. 

Acknowledgments 

We acknowledge Professor M. Cirillo for providing us with the photolithographic masks 
we used to fabricate the devices. Fruitful discussions with S. Pagano and C. Nappi, as well 
as the financial support of MURST COFINOO are also acknowledged. 



* Corresponding author; 
e-mail: ^iocar@sa. infn.it ; 
FAX: +3908965390 

1 N. Gr0nbech-Jensen, P. S. Lomdahl, and M. R. Samuelsen, Phys. Rev. B 43, 12799 (1991). 



15 



2 G. Carapella, Phys. Rev. B 63, 054515 (2001). 

^ G. Carapella and G. Costabile, Phys. Rev. Lett. 87, 077002 (2001). 

^ A. WallrafF, Y. Koval, M. Levitchev, M. V. Fistul, and A. V. Ustinov, J. Low. Temp. Phys. 
118, 543 (2000). 

^ S. Sakai, M. R. Samuelsen, and O. H. Olsen, Phys. Rev. B 36, 217 (1987). 

^ S. Pagano, C. Nappi, R. Cristiano, E. Esposito, L. Prunzio, L. Parlato, G. Peluso, G. Pepe, and 

U. Scotti Di Uccio, in Nonlinear Superconducting Devices and High Tc Materials, edited by R. 

D. Parmentier and N. F. Pedersen (World Scientific, Singapore, 1995). 
^ A. Benabdallah, J. G. Caputo, and A. C. Scott, Phys. Rev. B 54, 16139 (1996). 
8 A. Benabdallah, J. G. Caputo, and A. C. Scott, J. Appl. Phys. 88, 3527 (2000). 
^ E. Goldobin, A. Sterck, and D. Koelle, Phys. Rev. E 63, 031111 (2001). 

° V. P. Koshelets, P. N. Dmitriev, A. S. Sobolev, A. L. Pankratov, V. V. Khodos, V. L. Vaks, A. 

M. Baryshev, P. R. Wesselius, and J. Mygind, to be pubhshed in PHYSICA C (2001). 
^ T. Nagatsuma, K. Enpuku, K. Sucoka, K. Yoshida, and F. Irie, J. Appl. Phys. 58, 4412 (1985). 
^ V. P. Koshelets, S. V. Shitov, A. V. Shuchukin, L. V. Filippenko, J. Mygind, and A. V. Ustinov, 

Phys. Rev. B 56, 5572 (1997). 
^ A. Barone and G. Paterno, Physics and Applications of the Josephson Effect (John Wiley &z, 

Sons, New York, 1982). 

^ M. Cirillo, N. Gr0nbech-Jensen, M. R. Samuelsen, M. Salerno, and G. Verona Rinati, Phys. 

Rev. B 58, 12377 (1998). 
^ D. W. McLaughlin and A. C. Scott, Phys. Rev. A 18, 1652 (1978). 
^ A. Barone, J. Appl. Phys. 42, 2747 (1971). 

^ A. C. Scott and W. J. Johnson, Appl. Phys. Lett. 14, 316 (1969). 

^ S. Pace and U. Gambardella, J. Low. Temp. Phys. 62, 197 (1986). 

^ A. V. Ustinov, H. Kohlstedt, and P. Henne, Phys. Rev. Lett. 77, 3617 (1996). 

'° P. Cikmacs, M. Cirillo, V. Merlo, and R. Russo, IEEE Trans. Appl. Sup. 11, 99 (2001). 



16 



